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Abstract 

We study the optimal portfolio liquidation problem over a finite horizon in a limit 
order book with bid-ask spread and temporary market price impact penalizing speedy 
execution trades. We use a continuous-time modeling framework, but in contrast with 
previous related papers (see e.g. [24j and [25j). we do not assume continuous-time 
trading strategies. We consider instead real trading that occur in discrete-time, and 
this is formulated as an impulse control problem under a solvency constraint, including 
the lag variable tracking the time interval between trades. A first important result 
of our paper is to show that nearly optimal execution strategies in this context lead 
actually to a finite number of trading times, and this holds true without assuming ad hoc 
any fixed transaction fee. Next, we derive the dynamic programming quasi-variational 
inequality satisfied by the value function in the sense of constrained viscosity solutions. 
We also introduce a family of value functions converging to our value function, and 
which is characterized as the unique constrained viscosity solutions of an approximation 
of our dynamic programming equation. This convergence result is useful for numerical 
purpose, postponed in a further study. 
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1 Introduction 



Understanding trade execution strategies is a key issue for financial market practitioners, 
and has attracted a growing attention from the academic researchers. An important pro- 
blem faced by stock traders is how to liquidate large block orders of shares. This is a 
challenge due to the following dilemma. By trading quickly, the investor is subject to 
higher costs due to market impact reflecting the depth of the limit order book. Thus, 
to minimize price impact, it is generally beneficial to break up a large order into smaller 
blocks. However, more gradual trading over time results in higher risks since the asset 
value can vary more during the investment horizon in an uncertain environment. There 
has been recently a considerable interest in the literature on such liquidity effects, taking 
into account permanent and/or temporary price impact, and problems of this type were 
studied by Bertsimas and Lo [7], Almgren and Criss [T], Bank and Baum [5], Cetin, Jarrow 
and Protter [8], Obizhaeva and Wang [18], He and Mamayski [13], Schied an Schoneborn 
[25], Ly Vath, Mnif and Pham [17], Rogers and Singh [24j . and Cetin, Soner and Touzi [9], 
to mention some of them. 

There are essentially two popular formulation types for the optimal trading problem 
in the literature: discrete-time versus continuous-time. In the discrete-time formulation, 
we may distinguish papers considering that trading take place at fixed deterministic times 
(see [7]), at exogenous random discrete times given for example by the jumps of a Poisson 
process (see [22], [6]), or at discrete times decided optimally by the investor through an 
impulse control formulation (see [13] and [Uj). In this last case, one usually assumes the 
existence of a fixed transaction cost paid at each trading in order to ensure that strategies 
do not accumulate in time and occur really at discrete points in time (see e.g. [H] or |19j). 
The continuous-time trading formulation is not realistic in practice, but is commonly used 
(as in [8|, [25] or [2^), due to the tractability and powerful theory of the stochastic calculus 
typically illustrated by Ito's formula. In a perfectly liquid market without transaction cost 
and market impact, continuous-time trading is often justified by arguing that it is a limit 
approximation of discrete-time trading when the time step goes to zero. However, one may 
question the validity of such assertion in the presence of liquidity effects. 

In this paper, we propose a continuous-time framework taking into account the main 
liquidity features and risk/cost tradeoff of portfolio execution: there is a bid-ask spread 
in the limit order book, and temporary market price impact penalizing rapid execution 
trades. However, in contrast with previous related papers ( |25j or [24]), we do not as- 
sume continuous-time trading strategies. We consider instead real trading that take place 
in discrete-time, and without assuming ad hoc any fixed transaction cost, in accordance 
with the practitioner literature. Moreover, a key issue in line of the banking regulation 
and solvency constraints is to define in an economically meaningful way the portfolio value 
of a position in stock at any time, and this is addressed in our modelling. These issues 
are formulated conveniently through an impulse control problem including the lag variable 
tracking the time interval between trades. Thus, we combine the advantages of the stochas- 
tic calculus techniques, and the realistic modeling of portfolio liquidation. In this context, 
we study the optimal portfolio liquidation problem over a finite horizon: the investor seeks 
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to unwind an initial position in stock shares by maximizing his expected utihty from ter- 
minal liquidation wealth, and under a natural economic solvency constraint involving the 
liquidation value of a portfolio. 

A first important result of our paper is to show that that nearly optimal execution 
strategies in this modeling lead actually to a finite number of trading times. While most 
models dealing with trading strategies via an impulse control formulation assumed fixed 
transaction cost in order to justify a posteriori the discrete-nature of trading times, we 
prove here that discrete-time trading appear naturally as a consequence of liquidity fea- 
tures represented by temporary price impact and bid-ask spread. Next, we derive the 
dynamic programming quasi- variational inequality (QVI) satisfied by the value function in 
the sense of constrained viscosity solutions in order to handle state constraints. There are 
some technical difficulties related to the nonlinearity of the impulse transaction function 
induced by the market price impact, and the non smoothness of the solvency boundary. In 
particular, since we do not assume a fixed transaction fee, which precludes the existence of 
a strict supersolution to the QVI, we can not prove directly a comparison principle (hence 
a uniqueness result) for the QVI. We then consider two types of approximations by in- 
troducing families of value functions converging to our original value function, and which 
are characterized as unique constrained viscosity solutions to their dynamic programming 
equations. This convergence result is useful for numerical purpose, postponed in a further 
study. 

The plan of the paper is organized as follows. Section 2 presents the details of the model 
and formulates the liquidation problem. In Section 3, we show some interesting economical 
and mathematical properties of the model, in particular the finiteness of the number of 
trading strategies under illiquidity costs. Section 4 is devoted to the dynamic programming 
and viscosity properties of the value function to our impulse control problem. We propose 
in Section 5 an approximation of the original problem by considering small fixed tran- 
saction fee. Finally, Section 6 describes another approximation of the model with utility 
penalization by small cost. As a consequence, we obtain that our initial value function is 
characterized as the minimal constrained viscosity solution to its dynamic programming 
QVI. 

2 The model and liquidation problem 

We consider a financial market where an investor has to liquidate an initial position of 
y > shares of risky asset (or stock) by time T. He faces with the following risk/cost 
tradeoff: if he trades rapidly, this results in higher costs for quickly executed orders and 
market price impact; he can then split the order into several smaller blocks, but is then 
exposed to the risk of price depreciation during the trading horizon. These liquidity effects 
received recently a considerable interest starting with the papers by Bertsimas and Lo [7], 
and Almgren and Criss [l] in a discrete-time framework, and further investigated among 
others in Obizhaeva and Wang [18], Schied an Schoneborn [25], or Rogers and Singh [2^ 
in a continuous-time model. These papers assume continuous trading with instantaneous 
trading rate inducing price impact. In a continuous time market framework, we propose 
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here a more realistic modeling by considering that trading takes place at discrete points in 
time through an impulse control formulation, and with a temporary price impact depending 
on the time interval between trades, and including a bid-ask spread. 

We present the details of the model. Let (O, T ^ P) be a probability space equipped with 
a filtration F = {J^t)o<t<T satisfying the usual conditions, and supporting a one dimensional 
Brownian motion W on a finite horizon [0, T], T < oo. We denote by P = (P^) the market 
price process of the risky asset, by Xt the amount of money (or cash holdings), by Yj the 
number of shares in the stock held by the investor at time t, and by Qt the time interval 
between time t and the last trade before t. We set = (0, oo) and M^, = (—00, 0). 

• Trading strategies. We assume that the investor can only trade discretely on [0, T]. 
This is modelled through an impulse control strategy a = (r„, Cn)n>o- "^o < • • • < "^n • • • < 
are nondecreasing stopping times representing the trading times of the investor and Cn; 
n > 0, are ^t-„— measurable random variables valued in M and giving the number of stock 
purchased if Cn > or selled if Cn < at these times. We denote by A the set of trading 
strategies. The sequence {Tn,Cn) rnay be a priori finite or infinite. Notice also that we 
do not assume a priori that the sequence of trading times (?:„) is strictly increasing. We 
introduce the lag variable tracking the time interval between trades: 

Qt = inf {t-r„ : r„ < t}, t£ [0,r], 

which evolves according to 

@t = t-Tn, Tn<t<Tn+l, Qr„+, = 0, 71 > 0. (2.1) 

The dynamics of the number of shares invested in stock is given by: 

Yt = Yr„, Tn<t<Tn+l, Yr^^^=Y^- +Cn+1, n > 0. (2.2) 

ri + 1 

• Cost of illiquidity. The market price of the risky asset process follows a geometric 
Brownian motion: 

dPt = Ptibdt + adWt), (2.3) 

with constants b and a > 0. We do not consider a permanent price impact on the price, 
i.e. the lasting effect of large trader, but focus here on the effect of illiquidity, that is the 
price at which an investor will trade the asset. Suppose now that the investor decides at 
time t to make an order in stock shares of size e. If the current market price is p, and the 
time lag from the last order is 6, then the price he actually get for the order e is: 

Q{e,p,e) = pf{e,e), (2.4) 

where / is a temporary price impact function from R x [0,T] into U {oo}. We assume 
that the Borelian function / satisfies the following liquidity and transaction cost properties: 
(Hlf) /(0, 9) = 1, and /(., 9) is nondecreasing for all 9 G [0, T], 
(H2f) (i) /(e,0) = for e < 0, and (ii) /(e,0) = oo for e > 0, 
(H3f) Kb := sup(g g)gK._x[o,T] /(e,^') < 1 and := inf(e^0)giR*^x[o,T] f{e,9) > 1. 
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Condition (Hlf) means that no trade incurs no impact on the market price, i.e. Q{0,p, 9) 
= p, and a purchase (resp. a sale) of stock shares induces a cost (resp. gain) greater 
(resp. smaller) than the market price, which increases (resp. decreases) with the size of the 
order. In other words, we have Q{e,p,9) > (resp. <) p for e > (resp. <) 0, and Q{.,p,9) 
is nondecr easing. Condition (H2f) expresses the higher costs for immediacy in trading: 
indeed, the immediate market resiliency is limited, and the faster the investor wants to 
liquidate (resp. purchase) the asset, the deeper into the limit order book he will have to go, 
and lower (resp. higher) will be the price for the shares of the asset sold (resp. bought), with 
a zero (resp. infinite) limiting price for immediate block sale (resp. purchase). Condition 
(H2f) also prevents the investor to pass orders at consecutive immediate times, which is 
the case in practice. Instead of imposing a fixed arbitrary lag between orders, we shall see 
that condition (H2) implies that trading times are strictly increasing. Condition (H3f) 
captures a transaction cost effect: at time t, Pt is the market or mid-price, Ki,Pt is the bid 
price, KaPt is the ask price, and (kq — K,b)Pt is the bid-ask spead. We also assume some 
regularity conditions on the temporary price impact function: 

(Hcf) (i) / is continuous on R* x (0,T], 

df 

(ii) / is on Ml x [0, T] and x is bounded on Rt x [0, T]. 

A usual form (see e.g. |16j . |23j . [5]) of temporary price impact and transaction cost function 
/, suggested by empirical studies is 

f(e,9) = e^ltl''«g°(^)(K,l,>o + le=o + ^fele<o), (2.5) 

with the convention /(0,0) = 1. Here < k;, < 1 < k^, — ^fe is the bid-ask spread 
parameter, A > is the temporary price impact factor, and /3 > is the price impact 
exponent. In our illiquidity modelling, we focus on the cost of trading fast (that is the 
temporary price impact), and ignore as in Cetin, Jarrow and Protter and Rogers and 
Singh [231 the permanent price impact of a large trade. This last effect could be included 
in our model, by assuming a jump of the price process at the trading date, depending on 
the order size, see e.g. He and Mamayski [13j and Ly Vath, Mnif and Pham [17j . 

• Cash holdings. We assume a zero risk-free return, so that the bank account is constant 
between two trading times: 

Xt = Xr„, Tn<t<Tn+l, n > 0. (2.6) 

When a discrete trading AYt = Cn+i occurs at time t = r„+i, this results in a variation of 
the cash amount given by AXt := X^ — X^- = —/SYt.Q{/\Yt, Pt, G^- ) due to the illiquidity 
effects. In other words, we have 

= ~ Cn+l'3(Cn+l, -Pt„+i , 0^- ) 

n+1 n+1 

= X - - Cn+l-Pr„+i/(Cn+l>'^n+l - T"n), n > 0. (2.7) 

n + 1 

Notice that similarly as in the above cited papers dealing with continuous-time trading, 
we do not assume fixed transaction fees to be paid at each trading. They are practically 
insignificant with respect to the price impact and bid-ask spread. We can then not exclude 
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a priori trading strategies with immediate trading times, i.e. 6 - = r„_|_i — t„ = for 
some n. However, notice that under condition (H2f), an immediate sale does not increase 
the cash holdings, i.e. ^t„+i = X - = Xt-„, while an immediate purchase leads to a 
bankruptcy, i.e. -'^m+i = — oo. 

• Liquidation value and solvency constraint. A key issue in portfolio liquidation is 
to define in an economically meaningful way what is the portfolio value of a position on 
cash and stocks. In our framework, we impose a no-short sale constraint on the trading 
strategies, i.e. 

Yt > 0, < t < T, 

which is in line with the bank regulation following the financial crisis, and we consider the 
liquidation function L{x,y,p,9) representing the net wealth value that an investor with a 
cash amount x, would obtained by liquidating his stock position y > by a single block 
trade, when the market price is p and given the time lag 6 from the last trade. It is defined 
on R X M+ X R:; X [0, T] by 

L{x,y,p,e) = X + ypf{-y,9), 

and we impose the liquidation constraint on trading strategies: 

L{Xt,Yt,Pt,Qt) > 0, o<t<r. 

We have L{x,0,p,9) = x, and under condition (H2f)(ii), we notice that L(x,y,p, 0) = x 
for y > 0. We naturally introduce the liquidation solvency region: 

S = |(z,6') = (x,y,p,6l) G R X R+ X R; X [0,r] : y > and L(z,6') > o|. 

We denote its boundary and its closure by 

dS = dySUdiS and S = SUdS, 

where 

dyS = |(z,6') = (x,y,p,6') G R X R+ X R; X [0,T] 
diS = |(z,6i) = (x,y,p,6') G R X R+ X R; X [0,T] 
We also denote by Dq the corner line in dS: 

Do = {0} X {0} X R; X [0, T] = dyS n diS. 

• Admissible trading strategies. Given (t, z, 0) G [0, T] x 5, we say that the impulse 
control strategy a = (r„,,Cn)n>o is admissible, denoted by a G A{t,z.,9), \i = t — 9, r„ 
> t, n > 1, and the process {{Zs,Qs) = iXsYs,Ps,@s),t < s < T} solution to ([2J |) - ([2:2]) - 
<^M-UM-<^^i with an initial state (^^-,6^-) = {z,9) (and the convention that (^(,9*) 
= {z,9) if Ti > t), satisfies {Zs,@s) G [0,T] x S for all s G [t,T]. As usual, to alleviate 
notations, we omitted the dependence of {Z, Q) in (t, z, 9, a), when there is no ambiguity. 



y = and x = L{z,9) > o|, 
L{z,9) = 0y 
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theta=1 .5 



y : stock shares 



theta=0.5 



Do 



y : stock shares 



theta=1 



y : stock shares 
theta=0.1 





Do 
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y : stock shares 



Figure 1: Domain S in the nonhatched zone for fixed p = 1 and 9 evolving from 1.5 to 0.1. 
Here = 0.9 and f{e,9) = Kbexp(|) for e < 0. Notice that when 9 goes to 0, the domain 
converges to the open orthant M!^ x M^. 
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theta=1 




Figure 2: Lower bound of the domain S for fixed 6 = 1. Here Kh = 0.9 and f{e,6) = 
Kbexp(^) for e < 0. Notice that when p is fixed, we obtain the Figure 1. 



p=i 




Figure 3: Lower bound of the domain S for fixed p = 1 with f{e,6) = Kbexp(|) for e < 
and Kf, = 0.9. Notice that when 6 is fixed, we obtain the Figure 1. 
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Remark 2.1 Let {t,z,0) G [0, T] x S, and consider the impulse control strategy a = 
(Tm Cn)n>0) tq = t — 9 , Consisting in liquidating immediately all the stock shares, and then 
doing no transaction anymore, i.e. {ti,(i) = {t,—y), and = 0, n > 2. The associated 
state process (Z = {X, Y, P), Q) satisfies Xs = L{z, 9)^Ys = 0, which shows that L{Zs, Qg) 
= Xs = L{z, 9)>0,t<s<T, and thus a G A{t, z, 9) / 0. 



• Portfolio liquidation problem. We consider a utility function U from ]R_(_ into M, 
nondecreasing, concave, with C/(0) = 0, and s.t. there exists K > Q and 7 G [0, 1): 

(HU) < U{x) < Kx^, Vx G M+. 

The problem of optimal portfolio liquidation is formulated as 

v{t,z,9) = sup E[C/(Xt)], (t, z, 6') G [0, T] X cS, (2.8) 

where A£{t,z,9) = |a G A{t,z,9) : Yt = 0| is nonempty by Remark 12. li Notice 
that for a G Aii{t,z,9), Xt = L{Zt,Qt) > 0, so that the expectations in (j2.8|) . and the 
value function v are well-defined in [0, 00]. Moreover, by considering the particular strategy 
described in Remark 12.11 which leads to a final liquidation value Xt = L{z,9), we obtain 
a lower-bound for the value function; 

v{t,z,9) > U{L{z,9)), {t,z,9) e[0,T]xS. (2.9) 

Remark 2.2 We can shift the terminal liquidation constraint in Ai{t,z,9) to a terminal 
liquidation utility by considering the function Ul defined on S by: 

Ul{z,9) = U{L{z,9)), {z,9)eS. 

Then, problem (j2.8p is written equivalently in 



v{t,z,9) = sup E Ul{Zt,Qt) , {t,z,9) £ [0,T] X S. (2.10) 

aGA{t,z,e) L J 

Indeed, by observing that for all a G A£{t,z,9), we have E[C/(Xt)] = K[Ul{Zt,Qt)], and 
since Ai{t,z,9) C A{t,z,9), it is clear that v < v. Conversely, for any a G A{t,z,9) as- 
sociated to the state controlled process (Z, 0), consider the impulse control strategy a = 
a U (T, —Yt) consisting in liquidating all the stock shares Yt at time T. The correspond- 
ing state process {Z,Q) satisfies clearly: {Zs,@s) = {Zs,Qs) for t < s < T, and Xt = 
L{Zt,@t), Yt = 0, and so d G Ai{t,z,9). We deduce that EpLiZT^Qr)] = E[U{Xt)] 
< v{t,z,9), and so by arbitrariness of a in A{t,z,9), v{t,z,9) < v{t,z,9). This proves 
the equality v = v. Actually, the above arguments also show that sup^^^g^^^^^^^ U{Xt) = 
s^PaeA{t,z,e) Ul{Zt,Qt)- 

Remark 2.3 A continuous-time trading version of our illiquid market model with stock 
price P and temporary price impact / can be formulated as follows. The trading strategy 
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is given by a F-adated process r] = {r]t)o<t<T representing the instantaneous trading rate, 
which means that the dynamics of the cumulated number of stock shares Y is governed by: 

dYt = rjtdt. (2.11) 

The cash holdings X follows 

dXt = -i]tPtfivt)dt. (2.12) 

Notice that in a continuous-time trading formulation, the time interval between trades is 
G( = at any time t. Under condition (H2f), the liquidation value is then given at any 
time t by: 

L{Xt,Yt,Pt,0) = Xt, 0<t<T, 

and does not capture the position in stock shares, which is economically not relevant. On 
the contrary, by explicitly considering the time interval between trades in our discrete-time 
trading formulation, we take into account the position in stock. 

3 Properties of the model 

In this section, we show that the illiquid market model presented in the previous section 
displays some interesting and economically meaningful properties on the admissible trading 
strategies and the optimal performance, i.e. the value function. Let us consider the impulse 
transaction function T defined on M x M+ x M;; x [0, T] x M into R U {-co} x R x R;; by: 

r{z,9,e) = (x - epf{e,9),y + e,pj, 

for z = {x,y,p), and set r{z,9,e) = (r(z, ^, e), O) . This corresponds to the value of the 
state variable (Z, G) immediately after a trading at time t = r^+i of Cn+i shares of stock, 
i.e. (Zr^ , Gt-„ , i) = (r(Z - , G - , Cn+i), O) . We then define the set of admissible trans- 
actions: 

C(z,0) = |e G R : (r(z,^,e),0) G 5}, (z, e) G 5. 

This means that for any a = {TnXn)n>o S A{t,z,6) with associated state process {Z,Q), 
we have Cn G C(Z -, G^-), n > 1. We define the impulse operator TC by 

n(p{t,z,9) = sup ip{t,T{z,9,e),0), {t, z,9) e [0,T] x S. 
eeC{z,e) 

We also introduce the liquidation function of the (perfectly liquid) Merton model: 

Lm{z) = x + py, \/z = {x,y,p) eRxRxM.*^. 

For {t,z = {x,y,p),9) G [0,T] x S, we denote by (Z°'*'^, G°'*'^) the state process starting 
from (z, 9) at time t, and without any impulse control strategy: it is given by 

(^Z^Az^qOAB^ = (x,y,Pj^P,0 + s-t), t<s<T, 
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where P*'^ is the solution to (j2.3p starting fromp at time t. Notice that (Z^'*'^, B*''*'^) is the 
continuous part of the state process {Z,Q) controlled by a E A{t,z,9). The infinitesimal 
generator C associated to the process (Z*^'*'^, 0*^'*'^) is 

d^p dip 1 2 2d'^V> dcp 

We first prove a useful result on the set of admissible transactions. 

Lemma 3.1 Assume that (Hlf), (H2f) and (H3f) hold. Then, for all {z = {x,y,p),9) 
G S, the set C{z,6) is compact in M and satisfy 

C{z,e) C [-y,e{z,9)], (3.1) 

where —y < e{z, 0) < oo is given by 

f sup|eGM: ep/(e, 6*) < , if 9 > 
~ I 0, if9 = 0. 



e[z, 



For 9 = 0, (|3.ip becomes an equality : C{z,0) = [—y,0]. 

The set function C is continous for the Hausdorff metric, i.e. if {zn,9n) converges to 
(z, 9) in S , and (cn 

) is a sequence inC{zn,9n) converging to e, then e £ C{z,9). Moreover, 
i/ e G M I— > ef{e,9) is strictly increasing for 9 £ (0,T], then for (z = {x,y,p),9) E OlS 
with 9 > 0, we have e{z,9) = —y, i.e. C{z,9) = {—y}- 

Proof. By definition of the impulse transaction function T and the liquidation function L, 
we immediately see that the set of admissible transactions is written as 

C{z,9) = |e G R : X - ep/(e,6l) > 0, and y + e > o| 

= {eeR: epf{e,9) <x}n[-y,^) =: Ci{z,9)n[-y,oo). (3.2) 

It is clear that C{z, 9) is closed and bounded, thus a compact set. Under (Hlf) and (H3f), 
we have limg^oo ep/(e, 0) = oo. Hence we get e{z,9) < oo and Ci{z,9) C {—oo,e{z,9)]. 
From ([32]), we get dM]). Suppose 9 = 0. Under (H2f), using {z,9) £ S, we have Ci{z,9) 
= M_. From we get C{z,9) = [-y,0]. 

Let us now prove the continuity of the set of admissible transactions. Consider a 
sequence (zn = ixn,yn,Pn),9n) in S converging to {z,9) £ 5, and a sequence (en) in 
C{zn,9n) converging to e. Suppose first that 9 > 0. Then, for n large enough, 9n > 
and by observing that {z, 9, e) i-^ T{z, 9, e) is continuous on M x IR+ x M^j, x IR;!j_ x M, we 
immediately deduce that e G C{z, 9). In the case 9 = 0, writing Xn — en/(e„, 9n) > 0, using 
(H2f)(ii) and sending n to infinity, we see that e should necessarily be nonpositive. By 
writing also that yn + > 0, we get by sending n to infinity that y + e > 0, and therefore 
eeC(z,0) = [-y,0]. 

Suppose finally that e G M i— > ef{e,9) is increasing, and fix {z = {x,y,p),9) G OlS, 
with 9 > 0. Then, L{z,9) = 0, i.e. x = —ypf{—y,9). Set e = e{z,9). By writing that 
epf{e,9) < x = —ypf{—y,9), and e > —y, we deduce from the increasing monotonicity of 
e 1-^ epf{e, 9) that e = —y. □ 
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Remark 3.1 The previous Lemma implies in particular that C{z,0) C M_, which means 
that an admissible transaction after an immediate trading should be necessarily a sale. In 
other words, given a = {Tn,Cn)n>o £ -^(^i -2)^)1 (^,^,6) E [0, T] x S, if 0^- = 0, then < 
0. The continuity property of C ensures that the operator TC preserves the lower and upper- 
semicontinuity (see Appendix). This Lemma also asserts that, under the assumption of 
increasing monotonicity of e — > e/(e, 0), when the state is in the boundary L = 0, then the 
only admissible transaction is to liquidate all stock shares. This increasing monotonicity 
means that the amount traded is increasing with the size of the order. Such an assumption 
is satisfied in the example (j2.5p of temporary price impact function / for /5 = 2, but is not 
fulfilled for /5 = 1. In this case, the presence of illiquidity cost implies that it may be more 
advantageous to split the order size. 

We next state some useful bounds on the liquidation value associated to an admissible 
transaction. 

Lemma 3.2 Assume that (Hlf) holds. Then, we have for all {t,z,9) G [0,T] x S: 

< L{z,6) < Lm{z), (3.3) 
LM{T{z,e,e)) < Lm{z), VeGM, (3.4) 
sup L{Zs,es) < Lm(Z°'*'^), t<s<T. (3.5) 

a€A{t,z,e) 

Furthermore, under (H3f), we have for all [z = {x,y,p),6) G S, 

LM{^{z,e,e)) < Lm{z) -Tam.{Ka- 1,1 - iib)\e\p-, Ve G M. (3.6) 

Proof. Under (Hlf), we have f{e,6) < 1 for all e < 0, which shows clearly (13. 3p . Prom 
the definition of Lm and T, we see that for all e G M, 

LM(r{z,e,e))-LM{z) = ep(l-f{e,e)), (3.7) 

which yields the inequality (|3.4p . Fix some arbitrary a = (r^, Cn)n>o £ •^{t, z, 9) associated 
to the controlled state process (Z, 0). When a transaction occurs at time s = n > 1, 
the jump of Lm{Z) is nonpositive by (|3.4p : 

/^Lm{Zs) = LMiZr„)-LM{Z^~) = LM(r(Z-,e^-,Cn))-iA/(Z,-) < 0. 

We deduce that the process Lm{Z) is smaller than its continuous part equal to Lm{Z^'*'^), 
and we then get ()3.5p with ()3.3p . Finally, under the additional condition (H3f), we easily 
obtain inequality (j3.6p from relation (j3.7p . □ 

We now check that our liquidation problem is well-posed by stating a natural upper- 
bound on the optimal performance, namely that the value function in our illiquid market 
model is bounded by the usual Merton bound in a perfectly liquid market. 
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Proposition 3.1 Assume that (Hlf) and (HU) hold. Then, for all {t,z,9) G [0,r] x S, 
the family {Ul{Zt,Qt), ct E A{t,z,9)} is uniformly integrahle, and we have 



v{t,z,9) < vo{t,z) :-- 



E 



U L 



70,t,Z 



it,z,9)€[0,T]xS, 



where p is a positive constant s. t. 



P > 



7 



62 



1 - 7 2cj2 ■ 



(3.8) 



(3.9) 



Proof. From ()3.5p and the nondecreasing monotonicity of f7, we have for all (t, z, 9) € 
[0,T] X S: 

sup U{Xt) = sup Ul{Zt,Qt) < U{Lm{Z^''")), 
a&Ai{t,z,e) aeA(t,z,e) 

and all the assertions of the Proposition will follow once we prove the inequality (13. 8p . For 
this, consider the nonnegative function ip defined on [0, T] x S by: 

^{t,z,9) = e^^^~'^LM{zr = e''(^"*)(x+py)^ 

and notice that if is smooth on [0, T] x (5 \ Dq). We claim that for p > large enough, 
the function ip satisfies: 



dip dip 



dp > 0, on [0,r] X {S\Do). 



dt 89 

Indeed, a straightforward calculation shows that for all {t,z,9) G [0,T] x (S \ Dq): 
-^{t,z,9)-^{t,z,9)-Cipit,z,9) 



oPiT-t) 



Lm{z 



i7-2 



7 



Ap 



2 2 

y p 



(3.10) 



which is nonegative under condition (j3.9|) . 

Fix some (t, z, 9) G [0, T] x 5. If (z, 9) = (0, 0,p, 0) G Dqi then we clearly have t>o(i, z;, ^) 
= C/(0), and inequality (|3.8p is trivial. Otherwise, if {z,9) G 5 \ Dq, then the process 
(^o,i,^^0O,t,9) satisfy LA/(ZO'*'^ e^'*'^) > 0. Indeed, Denote by (Z*'^e*''') the process 
starting from (z, 0) at t and associated to the strategy consisting in liquidating all stock 
shares at t. Then we have {Zl'^ , Q^/) G 5 \ Dq for all s G [t, T] and hence Lni{Zl'^, 0s^)> 
for all s G [t,T]. Using ([M]) we get Lm{Zs'^'" > LM{Zi'\Q^/) > 0. 

We can then apply Ito's formula to (/^(s, Z^'*'^, 0s'*'^) between t and Tr = inf{s > 
t : > i?} AT: 



E[^(T«,Z°f,e?^*'^)] = ip{t,z)+E[p{^ + ^+Cip){s,Z'/'\er'^ 



^)ds 



< ip{t,z). 
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(The stochastic integral term vanishes in expectation since the integrand is bounded before 
Tr). By sending R to infinity, we get by Fatou's lemma and since ip(T, z, 6) = Lm{,z)^'- 



E 



We conclude with the growth condition (HU). □ 

As a direct consequence of the previous proposition, we obtain the continuity of the 
value function on the boundary dyS, i.e. when we start with no stock shares. 

Corollary 3.1 Assume that (Hlf) and (HU) hold. Then, the value function v is con- 
tinuous on [0, T] X dyS, and we have 

v{t,z,9) = U{x), yt£[0,T],{z,9) = {x,0,p,9) £dyS. 
In particular, we have v{t,z,0) = U{0) = 0, for all {t,z,9) £ [0, T] x Do- 
Proof. From the lower-bound (|2.9p and the upper-bound in Proposition 13.11 we have for 
ah {t, z, 9) £ [0, T] X S, 

U(x + ypf{-y,9)) <v{t,z,9) < E[C/(Lm(4'*''))] = ^U{x + yP!^^)]. 



These two inequalities imply the required result. □ 

The following result states the finiteness of the total number of shares and amount 
traded. 

Proposition 3.2 Assume that (Hlf) and (H3f) hold. Then, for any a = {Tn,Cn)n>o S 
A{t,z,9), {t,z,9) G [0,r] X S, we have 



^\Cn\ < OO, ^|Cn|-Pr„ < OO, and ^ ICnl^Vn/fCn,©^-) < OO, Q.S 



n>l n>l n>l 



Proof. Fix {t,z = {x,y,p),9) G [0,T] x S, and a = (r„,Cn)n>o G A(t,z,9). Observe 
first that the continuous part of the process Lm{Z) is Lm{Z^'^'^), and we denote its jump 
at time by ALa,/(.^t„) = L]\f{Zr,J — Lm{Z^-). From the estimates ()3.3p and (|3.6p in 
Lemma 13.21 we then have almost surely for all n > 1, 



< Lm(^.J = Lm(^^*'^) + 5^ALm(^.J 

k=l 
n 

< LM(^^*'^)-/^J]|GI^'r„ 



k=l 

where we set R = min(Ka — 1, 1 — k^) > 0. We deduce that for all n > 1, 

n 

y2\Ck\Pr, < - sup La/(Z°'*'^) = -{x + y sup P'.^P) < OO, a.s. 
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This shows the almost sure convergence of the series |Cn|-fT„- Moreover, since the price 
process P is continous and strictly positive, we also obtain the convergence of the series 
ICnl- Recalling that f{e,6) < 1 for all e < and 6 G [0,T], we have for all n > 1. 

n n n 

k=l k=l k=l 

n n 

< Y.^kPrj{Ck,&,^,)+2Y,\Ck\Pr,. (3.11) 
k=l k=l 

On the other hand, we have 

< LMiZrJ = Xr„+Yr„Pr„ 

n n 

= x-Y,CkPrJ{Ck,Qr-) + {y + Y.Ck)Pr^. 

k=l k=l 

Together with (|3.1ip . this implies that for all n > 1, 



T.\^k\Prj{Ck,e,-) < X + {y + J2\Ck\) sup P','P + 2Y,\Ck\Pr,. 

The convergence of the series \ Cn\PT„f {Cn, ©r") follows therefore from the convergence 
of the series En lCn| and Y,n \Cn\Pr„- ° 

As a consequence of the above results, we can now prove that in the optimal portfolio 
liquidation, it suffices to restrict to a finite number of trading times, which are strictly 
increasing. Given a trading strategy a = (T„,Cn)n>o S A, let us denote by N{a) the 
process counting the number of intervention times: 

Nt{a) = J^lr„<t, 0<t<r. 

n>l 

We denote by A'^{t, z, 9) the set of admissible trading strategies in Ai{t, z, 6) with a finite 
number of trading times, such that these trading times are strictly increasing, namely: 

A\{t,z,e) = ^a = {Tn,Cn)n>o^ Ai{t,z,e) : Nt {a) < oo, a.s. 

and Tn < Tn+i a.s., < n < NT{a) — l|. 
For any a = {Tn,Cn)n G A'}{t,z,6), the associated state process {Z,@) satisfies G - > 0, 

n + l 

i.e. (Z - ,G - ) eS* := \(z,e) G 5 : > ol. We also set OlS* = OlS n S*. 
Theorem 3.1 Assume that (Hlf), (H2f), (H3f), (Hcf) and (HU) hold. Then, we have 
v{t,z,0) = sup E[U{Xt)], {t,z,e) e [0,T] X S. (3.12) 

a£A'l{t,z,e) 

Moreover, we have 

v{t,z,e) = sup E[U{Xt)], {t,z,e) e[0,T]x {S\dLS), (3.13) 

aeA'i_^ {t,z,e) 

where A\ {t,z,e) = {q G A\{t,z,e) : (Z„e,) G {S\dLS),t< s < T}. 
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Proof. 1. Fix {t,z,9) G [0,T] x S, and denote by A'^{t,z,0) the set of admissible trading 
strategies in A£{t, z, 9) with a finite number of trading times: 

A\{t,z,9) = |a = (rfc, Cfc)fc>o G ^^(i, 0) : iVrla) is bounded a.s. |. 

Given an arbitrary a = {Tk,Ck)k>o G Ai{t,z,9) associated to the state process {Z,@) = 
{X, Y, P, G), let us consider the truncated trading strategy a^"^ = (r^, Cfe)fc<n U (t„+i, —Y - ), 

which consists in liquidating all stock shares at time Tn+i- This strategy a^"^ lies in 
Ai{t, z,9), and is associated to the state process denoted by e^")). We then have 

k>n+l 

Now, from Proposition 13.21 we have 

E Cfe-Prfc/(Cfc, ©T--) — ' ■^^^'^ n ^ oo. 

fc>n.+l 

Moreover, since < y - = 1V„ goes to 1^ = as n goes to infinity, by definition of ct E 
Ai{t, z,9), and recalling that / is smaller than 1 on IR_ x [0,T], we deduce that 




— > a.s. when n — > oo. 

This proves that 

x!p^ — > Xt a.s. when n — > oo. 

Prom Proposition EUl the sequence {U{X^^))n>i is uniformly integrable, and we can apply 
the dominated convergence theorem to get 

E[U{X^^)] — > E[U{Xt)], when n ^ oo. 

From the arbitrariness of a G A£{t, z,9), this shows that 

v{t, z, 9) < v\t, z,9):= sup E [U{Xt)] , 

a&A\{t,z,e) 

and actually the equality v = since the other inequality -u^ < f is trivial from the inclusion 
A^g{t,z,9) c Aiit,z,9). 

2. Denote by the value function in the r.h.s. of (|3.12|) . It is clear that < = v since 
A\{t,z,9) C A\{t,z,9). To prove the reverse inequality we need first to study the behavior 
of optimal strategies at time T. Introduce the set 

A\{t,z,9) = [a = {TkXk)k^A\{t,z,9) : : = T} < l}, 

and denote by u'' the associated value function. Then we have 'ifi < v^. Indeed, let a 
= {Tk,Ck)k be some arbitrary element in A''^{t,z,9), {t,z = {x,y,p),9) G [0,T] x S. If 
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a G A''^{t,z,6) then we have v^{t,z,6) > E Ul{Zt,@t) , where {Z,@) denotes the pro- 



cess associated to a. Suppose now that a i J^((t,z,e). Set 



maxjA; : < T}. 



Then define the stopping time r' := '^"^2'^ and the ^^-z -measurable random variable := 
argmax{e/(e, T -t^) : e> -Yr,^J. Define the strategy a' = (r^, Ck)k<m U (r', 1V„ - C) U 
(T, ^'). From the construction of a' , we easily check that a' E ^''(t, z, 6) and E C/l(2't, @t) 



< E 



Ul{Z^,Q'j,) where {Z',Q') denotes the process associated to a'. Hence, we get v'' > 



We now prove that > ifi . Let a = {Tk,Ck)k be some arbitrary element in A^^{t,z,9), 
(t, z = (x, y,p),9) G [0, T] xS. Denote by iV = Nxia) the a.s. finite number of trading times 
in a. We set m = inf{0 < A; < A^ — 1 : t^-^-i = r/j} and M = sup{m-|- 1 < k < N : = Tm} 
with the convention that inf = sup0 = A^-|-l. We then define a' = (r^, Cfc)o<fc<Ar-(A/-m)+i 
G ^by: 



i'^k^ C'k) 



inXk), 

{Tm = TM,J2k=mCk), 
{Tk 



M 



=m+l 



Ci) 



for < A; < m 
for k = m and m < N, 
for m+l<k<N-{M 
for /c = iV - (M - ?n) + 1 



m) and m < N, 



where r' = with f = max{Tfc : r^. < T}, and we denote by {Z' = {X' , Y', P), Q') the 
associated state process. It is clear that {Zg,Q'g) = {Zs,Qs) for t < s < r^, and so -^(Vy- 



(r')- 



e 



(r')- 



Moreover, since Tm = tm, we have - = for m + 1 < k < M. 



From Lemma l3.ll (or Remark 13. ip . this implies that Ck ^ for ?ti -|- 1 < < M, and so 
C7V-(A/-m)+i ~ Sfcf=m+i Cfe ^ 0. We also recall that immediate sales does not increase the 
cash holdings, so that X-^f. = X-^^ for m + \ < k < M . We then get 



Xrp 



> 



Xj' - 
Xj'. 



N-{M-m)+l 



Af-(A/-m)+l' ®(t')- 



-) 



Moreover, we have = y + J2k=i Ck = Y^ = 0. By construction, notice that Tq < . . . < 
r^_l_]^. Given an arbitrary a G A^{t, z,d), we can then construct by induction a trading 
strategy a' G A^{t, z, 9) such that X'j, > Xt a.s. By the nondecreasing monotonicity of the 
utility function [/, this yields 

E[C/(Xt)] < nU{X'T)] < v'{t,z,9), 

and we conclude from the arbitrariness of a G A'^{t, z, 9): < v^, and thus v = = = 



3. Fix now an element (t, z, 9) G [0, T] x (5 \ OlS), and denote by v+ the r.h.s of (I3.13|) . It 
is clear that v > v+. Conversely, take some arbitrary a = {rk, Ck)k G A^^{t, z, 9), associated 
with the state process {Z, 0), and denote by = NT{a) the finite number of trading times 
in a. Consider the first time before T when the liquidation value reaches zero, i.e. = 
mi{t < s < T : L{Zs, ©s) = 0} A T with the convention inf = oo. We claim that there 
exists 1 < m < iV + 1 (depending on uj and q) such that = Tm., with the convention 
that m = -|- 1, tn+i = T if = T. On the contrary, there would exist 1 < A: < such 
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that Tk < < Tfc+i, and L{Zt-'^,@j 



0. Between and Tk+i, there is no trading, and 



so {Xs, Ys) = (Xt-^,, IVJ, = s - Tfc for < s < Tfc+i. We then get 



^Tfc + ^TkPsfi -Yr^,S - Tk), Tk<S< Tfc+l. 



(3.14) 



Moreover, since < L(Zt-^, Qt^) = Xr^, and L{ZT-a, ©t") = 0, we see with (|3.14p for s = r" 
that Yrf. Pt<^ / ( — IVfc 1 t"^ — Tk) should necessarily be strictly negative: 1^^ Pt-^ f{ — Yr^. ,T°' — Tk) 
< 0, a contradiction with the admissibility conditions and the nonnegative property of /. 

We then have = for some 1 < ?n < + 1. Observe that if m < A^, i.e. 
L{Zt-^,Qt-^) = 0, then U{L{Zt,Qt)) = 0. Indeed, suppose that Yr^ > and m < N. 
From the admissibility condition, and by Ito's formula to L{Z, Q) in (j3.14p between r" and 
T~+i, we get 



< L(z - ,9^- ; 

'm+1 'fe+1 



Y P 



Li^Zf'^ , ©T-a ) 

/3(y,„, s - T^)ds + af{-Yr,,s- T„,)dWs\ , (3.15) 



where P(y,e) = bf(-y,e) + ^(-y,0) is bounded on M+ x [0,T] by (Hcf)(ii). Since the 
integrand in the above stochastic integral w.r.t Brownian motion W is strictly positive, 
thus nonzero, we must have r" = Tm+i- Otherwise, there is a nonzero probability that the 
r.h.s. of (I3.15P becomes strictly negative, a contradiction with the inequality (I3.15p . 

Hence we get = 0, and thus L[Z - , © - ) = X^-^ = 0. From the Markov feature 

m + l m + 1 

of the model and Corollarv 13.11 we then have 



E 



C/fL(ZT,©T) 



< viTm,Zr^,QrJ = U {X^J 



0. 



Since U is nonnegative, this implies that U[L{Zt,@t)) = 0. Let us next consider the 
trading strategy a' = (t^, Cfc)o<fe<(m-i) ^ consisting in following a until time r", and 
liquidating all stock shares at time = Tm~i, and defined by: 



(.Tk^ C'k, 



{Tk, Ck 
Tm-1, 



for < k < m 
for k = m — 1, 



and we denote by (Z',©') the associated state process. It is clear that (Zg,©^) = {Zs,Qs) 
for t < s < Tm-i, and so L(Z^,©'g) = L{Zs,Qs) > for i < s < Tm~i- The liquidation 
at time r^-i (for m < N) yields Xr^_-^ = L{Z - , © - ) > 0, and l^r„i_i = 0. Since 

m — 1 m — 1 

there is no more trading after time Tm-i, the liquidation value for r^-i < s < T is given 
by: L{Zs,Qs) = Xr„^_^ > 0. This shows that a' £ ^^^(t,z,0). When m = + 1, we 
have a = a', and so X^ = L{Zlp, ©^) = L{Zt, ©r) = Xt- For m< N,we have U{X^) = 
U{L{Z^,e'T)) > = [/(L(Zt,0t)) = U{Xt). We then get U{X^) > U{Xt) a.s., and so 

E[U{Xt)] < E[U{X^)] < v+{t,z,9). 



We conclude from the arbitrariness of a G ■A'^it, z, 9): v < and thus v = v+. 



□ 
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Remark 3.2 If we suppose that the function e G M i-^ ef{e,6) is increasing for 9 G (0,T], 
we get the value of v on the bound 5^5*: v{t,z,9) = C/(0) = for {t,z = {x,y,p),9) S 
[0, T] X BlS* . Indeed, fix some point {t,z = {x,y,p),6) £ [0,T] x OlS*, and consider 
an arbitrary a = {Tk,Ck)k £ A^^{t,z,6) with state process (Z,©), and denote by N the 
number of trading times. We distinguish two cases: (i) If ri = t, then by Lemma |3. 11 the 
transaction is equal to —y, which leads to = 0, and a liquidation value L{Zt-^, Qti) 
= = L{z,6) = 0. At the next trading date T2 (if it exists), we get = y_- = 
with liquidation value L{Z^- ,Q^-) = 0, and by using again Lemma l3.ll we see that after 
the transaction at T2, we shall also obtain Xt-2 = = 0. By induction, this leads at the 
final trading time to X^j^ = Yr^^ = 0, and finally to Xt = Yt = 0. (ii) If ri > t, we claim 
that y = 0. On the contrary, by arguing similarly as in (|3.15p between t and , we have 
then proved that any admissible trading strategy a G A'^{t, z,9) provides a final liquidation 
value Xt = 0, and so 

v{t,z,9) = U{0) = 0, y{t,z,9)€[0,T]xdLS*. (3.16) 

Remark 3.3 The representation (j3.12p of the optimal portfolio liquidation reveals inte- 
resting economical and mathematical features. It shows that the liquidation problem in 
a continuous-time illiquid market model with discrete-time orders and temporary price 
impact with the presence of a bid-ask spread as considered in this paper, leads to nearly 
optimal trading strategies with a finite number of orders and with strictly increasing trading 
times. While most models dealing with trading strategies via an impulse control formulation 
assumed fixed transaction fees in order to justify the discrete nature of trading times, 
we prove in this paper that discrete-time trading appears naturally as a consequence of 
temporary price impact and bid-ask spread. 

The representation (|3.13p shows that when we are in an initial state with strictly posi- 
tive liquidation value, then we can restrict in the optimal portfolio liquidation problem to 
admissible trading strategies with strictly positive liquidation value up to time The 
relation (13.160 means that when the initial state has a zero liquidation value, which is not 
a result of an immediate trading time, then the liquidation value will stay at zero until the 
final horizon. 

4 Dynamic programming and viscosity properties 

In the sequel, the conditions (Hlf), (H2f), (H3f), (Hcf) and (HU) stand in force, and 
are not recalled in the statement of Theorems and Propositions. 

We use a dynamic programming approach to derive the equation satisfied by the value 
function of our optimal portfolio liquidation problem. Dynamic programming principle 
(DPP) for impulse controls was frequently used starting from the works by Bensoussan 
and Lions [3], and then considered e.g. in [28], [20], [T7] or |26j . In our context (recall the 
expression (|2.10p of the value function), this is formulated as: 
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Dynamic programming principle (DPP). For all {t,z,0) £ [0,T] x S, we have 

v{t,Z,e) = sup E[v{T,Zr,er)], (4.1) 

where r = r(a) is any stopping time valued in [t, T] eventually depending on the strategy 
a in (j4.1|) . More precisely we have : 

(i) for all a G A{t, z, 9), for all r G Tt^T, the set of stopping times valued in [t, T]: 

E[v{T,Zr,er)] < v{t,Z,6) (4.2) 

(ii) for all e > 0, there exists £ A{t, z, 6) s.t. for all r G Tt^x- 

vit,z,e)-e < EKT,z^,e^;)], (4.3) 

with (Z^,0^) the state process controlled by a^. 

The corresponding dynamic programming Hamilton-Jacobi-Bellman (HJB) equation is 
a quasi-variational inequality (QVI) written as: 



mm 



dv dv 



0, in[0,r)x5, (4.4) 



together with the relaxed terminal condition: 

min [v-Ul , v-Hv] = 0, in {T} x S. (4.5) 

The rigorous derivation of the HJB equation satisfied by the value function from the 
dynamic programming principle is achieved by means of the notion of viscosity solutions, 
and is by now rather classical in the modern approach of stochastic control (see e.g. the 
books [llj and [21]). There are some specificities here related to the impulse control and 
the liquidation state constraint, and we recall in Appendix, definitions of (discontinuous) 
constrained viscosity solutions for parabolic QVIs. The main result of this section is stated 
as follows. 

Theorem 4.1 The value function v is a constrained viscosity solution to ()4.4p -( l43]) . 

Proof. The proof of the viscosity supersolution property on [0, T) x S and the viscosity 
subsolution property on [0,T) x S follows the same lines of arguments as in ^7\, and is 
then omitted here. We focus on the terminal condition (14. 5p . 

We first check the viscosity supersolution property on {T} x S. Fix some {z, 9) G 5, 
and consider some sequence {tk, Zk,9k)k>i in [0, T) x S, converging to {T,z,9) and such 
that limfc v{tk, Zk,9k) = f*(T', z, 9). By taking the no impulse control strategy on [tfc, T], we 
have 

v{tk,zk,9k) > E[c/i(4'*'=-~^e^*'=•''=)]. 
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bince converges a.s. to {z, 9) when k goes to infinity by continuity of 

^^o,i,2^ 0O,t,6»^ in its initial condition, we deduce by Fatou's lemma that 

v,{T,z,e) > UL{z,e). (4.6) 

On the other hand, we know from the dynamic programming QVI that v > 7iv on [0, T) x 5, 
and thus 

v{tk,zk,6k) > Hv{tk,zk,dk) > TCv^{tk,Zk,Ok), V/c > 1. 
Recahing that Wv* is Isc, we obtain by sending k to infnity: 

v4T,z,e) > nv4T,z,e). 

Together with (|4.6|) . this proves the required viscosity supersolution property of (|4.5p . 

We now prove the viscosity subsolution property on {T} x S, and argue by contradiction 
by assuming that there exists (z, 9) £ S such that 

mm[v*{T,z,e) -UL{z,e) , v*{T,z,e) -nv*{T,z,e)] := 2e > 0. (4.7) 

One can find a sequence of smooth functions {f"')n>o on [0, T] x S such that if"' converges 
pointwisely to v* on [0, T] x 5 as n — > oo. Moreover, by (|4.7|) and recalhng that Tiv* is use, 
we may assume that the inequality 

min [99" -Ul , - Hv*] > e, (4.8) 

holds on some bounded neighborhood i?" of (T, z, 9) in [0, T] x 5, for n large enough. Let 
(tfc, Zk-, Ok)k>i be a sequence in [0, T)xS converging to (T, z, 6) and such that lim/j ^;(tfc, -Zfci ^fc) 
= v*{T, z, ^). There exists 5" > such that := [tfc, T] x B{zk,S'') x (^(^^ - ,5", ^fc + 5") n 

[0, T]^ C -B" for all A; large enough, so that (|4.8p holds on i?^. Since f is locally bounded, 
there exists some 77 > such that \v*\ < 77 on i?'". We can then assume that 99" > —2r] on 
B^. Let us define the smooth function (p'^ on [0, T) x 5 by 



(^"(t,z,e)+4?7 



1^ - Zk\ 
■n\2 



and observe that 



{v*-rk)it,z,0) < -7?, 



(5",0fc + (5'^)n [0,T] . Since 



5t 



(4.9) 



-00 as 



C^l > 0, on S^. 



(4.10) 



for (t,z,^) G [tfc,r] xaS(zfc,5") X 
t — > r, we have for A: large enough 

_ ^ _ ^ 

dt de 

Let a'^ = {'T'j,Cj)j>i be a ^—optimal control for v{tk, Zk,Ok) with corresponding state 
process {Z^, e''), and denote by cr^ = inf{s > tk : (^^ Oj) ^ 5^} A rf A T. From the 
DPP (USD, this means that 

1 



v{tk,Zk,9k) 



k 



< E 



+ E 



+ E 



^;(rf,r(Z(';,)_,ef^,)_,Cf),0) 



(4.11) 
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Now, by applying Ito's Lemma to ifni^i ©s) between tk and a^, we get from (j4.8p - (j4.9 

dnni), 



> E 



'-a^<{TfAT) 



+ E 



k 



1,.<,. (v* (r^ r(Z(^^,). , , ), 0) + e 



Together with (j4.1ip . this imphes 

(fl{tk,Zk,Ok) > v{tk,Zk,Ok) - ^ + e At]. 

Sending k, and then n to infinity, we get the required contradiction: v*(T, z, 9) > v*{T, z, 6)+ 
£ Arj. □ 

Remark 4.1 In order to have a complete characterization of the value function through 
its HJB equation, we need a uniqueness result, thus a comparison principle for the QVI 
()4.4p - (|4.5p . A key argument originally due to Ishii [T3] for getting a uniqueness result for 
variational inequalities with impulse parts, is to produce a strict viscosity supersolution. 
However, in our model, this is not possible. Indeed, suppose we can find a strict viscosity 
Isc supersolution w to ()4.4p . so that {w — Ti.w)(t, z,6) > on [0, T) x S. But for z = 
{x, y,p) and 6 = 0, we have T{z, 0, e) = {x,y + e,p) for any e C{z, 0). Since G C{z, 0) we 
have 'Hw{t, z^Q) = supgg[_j^ q] x, y + e,p, 0) > w{t,z,Qi) > Tiwit, z,0), a contradiction. 
Actually, the main reason why one cannot obtain a strict supersolution is the absence of 
fixed cost in the impulse function T or in the objective functional. 



5 An approximating problem with fixed transaction fee 

In this section, we consider a small variation of our original model by adding a fixed 
transaction fee e > at each trading. This means that given a trading strategy a = 
iTn,Cn)n>o, the Controlled state process {Z = {X,Y, P),@) jumps now at time Tn+i, by: 



(^r„+,,e,„^j = r,(z ,e ,Cn+i),o , (5.1) 

\ 'n + l 'n + 1 / 

where is the function defined on M x M+ x M!; x [0, T] x E into MU {-cx)} x M x by: 
T,{z,e,e) = T{z, 9, e)-{e, 0,0) = (x-epf{e,e)-e,y + e,p), 



for z = {x, y,p). The dynamics of {Z, G) between trading dates is given as before. We also 
introduce a modified liquidation function defined by: 



Le{z,l 



max[x, L{z,e)-e], {z,e) = {x,y,p,0) £R xR+ xR* x [0,T]. 



The interpretation of this modified liquidation function is the following. Due to the presence 
of the transaction fee at each trading, it may be advantageous for the investor not to 
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liquidate his position in stock shares (which would give him L{z, 9) — e), and rather bin his 
stock shares, by keeping only his cash amount (which would give him x). Hence, the investor 
chooses the best of these two possibilities, which induces a liquidation value Lir{z,6). 

We then introduce the corresponding solvency region with its closure Sg = Ss D dSg, 
and boundary dSs = dySe U BlSs'- 



Se = 




dySe = 




OlSs = 





G M X M , X Ml X 



: L,{z,e)=Oy 



We also introduce the corner lines of dS^- For simplicity of presentation, we consider a 
temporary price impact function / in the form: 

fie,0) = f(^) = exp (A^ ) (Kale>0 + le=0 + K6le<0 )10>O, 



where < k;, < 1 < k^, and A > 0. A straightforward analysis of the function L shows 
that y 1-^ L{x,y,p,6) is increasing on [0,0/A], decreasing on [6/X,oo) with L{x,0,p,6) = 
X = L{x, oo,p,6), and maXy>o L{x, y,p, 9) = L{x,6 / X^p^ 0) = x + p^f{—l/\). We first get 
the form of the sets C{z, 6): 

e{z,e) = [-y,e{z,e)], 

where the function e is defined in Lemma 13.11 We then distinguish two cases: (i) If 
p^/(— 1/A) < e, then Lg{x,y,p,9) = x. (ii) If p^/(— 1/A) > e, then there exists an unique 
yi{p,9) G (0,e/A] and y2{p,0) £ [e/X,oo) such that L{x,yi{p,9),p,e) = L{x,y2ip,9),p,9) 
= X, and Le{x,y,p,9) = x for y e [0,yi{p,9))U{y2{p,9),oo), Le{x,y,p,9) = L{x,y,p,9)-e 
for y G [yi{p,9),y2{p,9)]. We then denote by 

Do = {0} X {0} X m; X [0, T] = dySe n OlSs, 
Di,e = {{0,yi{p,9),p,9):pjf{=^)>e, 9e[0,T]Y 

D2,s = {{0,y2{p,9),p,9):pjf{=^)>e, 9e[0,T]]. 

Notice that the inner normal vectors at the corner lines Di ^ and 1)2, e form an acute angle 
(positive scalar product), while we have a right angle at the corner Dq. 

Next, we define the set of admissible trading strategies as follows. Given {t,z,9) G 
[0, T] X cSj, we say that the impulse control a is admissible, denoted by a G A^(t,z,9), if 
tq = t — 9, Tn > t, n > 1, and the controlled state process {Z^,Q) solution to (|2.ip - (|2.2p - 
(I23])-(IM1)-(I5TI]), with an initial state {Z^^,Qt-) = {z,9) (and the convention that {Zf^Qt) 
= {z,9) if Ti > t), satisfies {Zg,Qs) G [Oi^] ^ for all s G [t,T]. Here, we stress the 
dependence of Z^ = {X^ , Y, P) in e appearing in the transaction function , and we notice 
that it affects only the cash component. Notice that A'^{t, z, 9) is nonempty for any (t, z, 9) 
G [0,T] X Indeed, for (z = {x,y,p),9) G S^, i.e. Lg{z,9) = max{x, L{z, 9) — e) > 0, 
we distinguish two cases: (i) if x > 0, then by doing none transaction, the associated state 
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epsilon=0.1 



epsilon=0.2 



D 




y : stock shares 
epsilon=0.3 



2,£ 




y : stock shares 



2,E 

\ 




y D ^ 




y : stock shares 
epsilon=0.4 



/ 



y : stock shares 



Figure 4: Domain in the nonhatched zone for fixed p = 1 and 6=1 and e evolving from 
0.1 to 0.4. Here Kb = 0.9 and /(e, 9) = K^exp (|) for e < 0. Notice that for e large enough, 
5^ is equal to open orthant Ml x Ml. 
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theta=1 




Figure 5: Lower bound of the domain for fixed 9 = 1 and /(e, 9) = K^exp (|) for e < 0. 
Notice that when p is fixed, we obtain the Figure 4. 



p=i 




Figure 6: Lower bound of the domain for fixed p = 1 and e = 0.2. Here Kh = 0.9 and 
f{e,6) = Kftexp (I) for e < 0. Notice that when 6 is fixed, we obtain the Figure 4. 
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process {Z^ = {X^,Y, P), Q) satisfies = x>0,t<s<T, and thus this zero transaction 
is admissible; (ii) if L{z,9) — e > 0, then by hquidating immediately all the stock shares, 
and doing nothing more after, the associated state process satisfies = L{z, 6) — e, Yg = 
0, and thus L^{Zg, Qg) = > 0, t < s < T, which shows that this immediate transaction 
is admissible. 

Given the utility function U on M_^, and the liquidation utility function defined on 
by Ul^{z,9) = U{Lir{z,9)), we then consider the associated optimal portfolio liquidation 
problem defined via its value function by: 

v,{t,z,e) = sup E[UlAZ^,Qt)], it,z,6) e[0,T]xSs. (5.2) 

Notice that when e = 0, the above problem reduces to the optimal portfolio liquidation 
problem described in Section 2, and in particular vo = v. The main purpose of this section 
is to provide a unique PDE characterization of the value functions v^, £ > 0, and to prove 
that the sequence {vs)e converges to the original value function u as e goes to zero. 

We define the set of admissible transactions in the model with fixed transaction fee by: 
C,{z,9) = [eeR: (T,{z,e,e),o) eSe}, {z,e)eSe. 
A similar calculation as in Lemma |3. II shows that for (z = {x,y,p),9) G S^, 



Ceiz,9) 



-y,e^{z,9)], if ^ > or X > e, 
0, if 6* = and x < e, 



where e{z, 9) = sup{e G M : epf{e/9) < x — e} ii 9 > and e{z, 0) = if x > e. Here, the 
set [—y,e£{z,9)] should be viewed as empty when e{z,9) < y, i.e. x + pyf{—y/9) — e < 0. 
We also easily check that Cg is continuous for the Hausdorff metric. We then consider the 
impulse operator Tie by 

new{t,z,9) = sup w{t,re{z,9,e),0), {t,z,9) e[0,T]xSe, 

for any locally bounded function w on [0, T] x 5e, with the convention that 7iew{t, z, 9) = 
— oo when Ce(z, 9) = 0. 

Next, consider again the Merton liquidation function and observe similarly as in 
(fHTD that 

LM{Te{z,9,e))-LM{z) = ep(l - f {e,9)) - e 

< -e, \/{z,9) gSs, e e M. (5.3) 

This implies in particular that 

HsLm < Lm on Se. (5.4) 

Since < Lm, we observe from (jS.Sp that if {z,9) £ A4 := {{z,9) £ Sg : Lm{z) < e}, 
then Ce(z, 9) = 0. Moreover, we deduce from ()5.3p that for all a = {tu, Cn)n>o £ -^^{i-, z, 9) 
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associated to the state process {Z,Q), {t,z,9) £ [0, T] x 5^: 

n>0 

< LM{Z^''')-eNT{a), 

where we recall that Nxia) is the number of trading times over the whole horizon T. This 
shows that 

Nxia) < -Lm{z!^'^'^) < oo a.s. 

In other words, we see that, under the presence of fixed transaction fee, the number of 
intervention times over a finite interval for an admissible trading strategy is finite almost 
surely. 



The dynamic programming equation associated to the control problem (j5.2p is 

dw dw 



mm 



Cw , w — HeW 



0, in[0,T)xcS„ (5.5) 



dt do 

min \_w — Ul^ , w — TLsw\ = 0, in {T} x Se- (5-6) 

The main result of this section is stated as follows. 

Theorem 5.1 (1) The sequence {ve)e is nonincreasing, and converges pointwise on [0,T] x 
{S \ diS) towards v as e goes to zero. 

(2) For any e > 0, the value function i;^ is continuous on [0,T) x S^, and is the unique (in 
[0,r) X Se) constrained viscosity solution to ()5.5p - (|5.6p . satisfying the growth condition: 

\vs{t,z,d)\ < K{l + LM{zr), y{t,z,d)e[0,T]xSs, (5.7) 

for some positive constant K , and the boundary condition: 

lim Ve{t\z',e') = v(t,z,9) 
{t',z',e')^it,z,e) 

= f/(0), y{t,z = {0,0,p),e)e[0,T]xDo. (5.8) 



We first prove the convergence of the sequence of value functions (ve)- 
Proof of Theorem (1). 

Notice that for any < ei < £2, we have L^^ < L^^ < L, A^^{t,z,9) C A'^^{t,z,0) C 
A{t,z,9), for t e [0,T], {z,e) G C Se, C 5, and for a £ A''^{t,z,e), Le^{Z^\Q) < 
Le2iZ^^ ,Q) < Le,{Z^^ ,@) < L{Z,Q). This shows that the sequence (ve) is nonincreasing, 
and is upper-bounded by the value function v without transaction fee, so that 

limveit,z,0) < v{t,z,9), V(t, z, 6*) G [0, T] x 5. (5.9) 

Fix now some point {t,z,6) £ [0, T] x (5 \ OlS). From the representation (j3.13p of 
v{t, z,9), there exists for any n > 1, an l/?i-optimal control = (t;["\ ci"^)fc G {t, z, 9) 
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with associated state process = (X("),y("),P),e(")) and number of trading times 

E\U{X^^)] > v{t,z,e)--. (5.10) 

We denote by (Z-^'W, G^")) = P), G^")) the state process controlled by a(") in 

the model with transaction fee e (only the cash component is affected by e), and we observe 
that for all t < s < T, 

Xf'W = Xi")-eA^i") / as e goes to zero. (5.11) 
Given n, we consider the family of stopping times: 

4") = inf {s > t : L(Zf'("), ei")) < e} A T, e > 0. 
Let us prove that 

limcri") = T a.s. (5.12) 

Observe that for < ei < £2, '^"^ < and so L(Zf '^"\g,) < L(Zf '^"\g,) 

for t < s < T. This implies clearly that the sequence {ai^^)s is nonincreasing. Since this 

(n) (n) 

sequence is bounded by T, it admits a limit, denoted by cJq = linigio T cr^ . Now, by 
definition of ai"'\ we have L{Z^f^\ G^"\) < e, for all e > 0. By sending s to zero, we then 
get with dSH]): 



^0 '^0 



Recalling the definition of (t, ^, 6), this implies that Cq = r^" for some /c G {1, . . . , A^^'"^ + 



1} with the convention t^(L , = T. If /c < A^("), arguing as in (j3.15p . we get a contradiction 



(n) 

with the solvency constraints. Hence we get = T. 



Consider now the trading strategy a^'^") g A consisting in following a^") until time cie"^ 
and liquidating all the stock shares at time al , i.e. 

^ (,W,cW)l^^^^^^„,U(aH,-i;(.),-). 

We denote by (Z^'{") = (X^-C")^ y^-C")^ p)^ 0e.("-)) the associated state process in the market 
with transaction fee e. By construction, we have for all t < s < cJe"'': L(Zs'^"\ G^'^"^) = 
L(Zs'^"\ gI"^) > e, and thus L£(Zf G^'^"^) > 0. At the transaction time cj^J^\ we then 
have = L{ztl , ^tl ) " ^ = ^(^^l),- ' t),- ) " = After time 4"^ , 

there is no more transaction in a'''("\ and so 

If'W = = > 0, (5.13) 

y^.,(n) ^ ^.,(n) ^ 0^ 4")<S<r, (5.14) 
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and thus ^^(Zl'^"^ e^'^"^) = Xf'^"^ > for a^^ < s <T. This shows that q^'(") hes in 
A^{t, z, 6), and thus by definition of v^: 



ve{t,z) > E[C/z.,(Z^'('^\e^("))]. (5.15) 



Let us check that given n, 

limL,(Z^'("\e^(")) = X!p\ a.s. (5.16) 

To alleviate notations, we set N = N!j?^ the total number of trading times of q^"). If the 
last trading time of a^"^ occurs strictly before T, then we do not trade anymore until the 
final horizon T, and so 

4") = XW, and y^") = = 0, on {r^v < T}. (5.17) 

By (|5.12|) . we have for e smah enough: a^"'^ > tn, and so X'^'!"^} = X^'jf'\ Y^lj^\ = Yr^^^ 
= 0. The final liquidation at time ai'^^ yields: 4'^"^ = X'^j'^) = X'^}^}_ — e = Xr'lf''^ — e, 
and y^'^"^ = Y'p'^ = 0. We then obtain 

L,(4'("\ e^(")) = max (4("\ L(Z^'("\ G^(")) - e) 

= 4(") = X^;(")-e on {r^ <T} 
= X^") - (1 + N)e on {tn < T}, 

by (j5.1ip and (|5.17p . which shows that the convergence in (j5.16p holds on {tn < T}. If the 
last trading of a^'"^ occurs at time T, this means that we liquidate all stock shares at T, 
and so 

x^") = L(z^"J,e$?2), y^") = o on {tn = t}. (5. is) 

On the other hand, by (|5.13p - (j5.14p . we have 

L,(4'("\e^(")) = x^'^'^) = L(z(:)„,_,e^';i,_)-6 

— > L{z!^},Q^^l) as e goes to zero, 

by (|5.12p . Together with (j5.18p . this implies that the convergence in (j5.16p also holds on 
{ttv = r}, and thus almost surely. Since < < L, we immediately see by Proposition 
I3.1l that the sequence {Ul^ (^Z^^^\ G)^''"'*) , e > 0} is uniformly integrable, so that by sending 
e to zero in (|5.15p and using (|5.16p . we get 

limve{t,z,e) > E[f/(4^)] >v{t,z)-^, 

from ()5.10p . By sending n to infinity, and recalling (15. 9p . this completes the proof of 
assertion (1) in Theorem 15. 1[ □ 
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We now turn to the viscosity characterization of v^- The viscosity property of is 
proved similarly as for and is then omitted. From Proposition 13.11 and since < f ^ 

< u, we know that the value functions Ve lie in the set of functions satisfying the growth 
condition in (j5.7p . i.e. 

g,([0,T]x5,) = {«;:[0,r]x5,-M, sup 1""})''^^ < oo|. 

The boundary property (|5.8p is immediate. Indeed, fix {t,z = {x,0,p),9) G [0,T] x dyS^, 
and consider an arbitrary sequence {tn,Zn = {xn,yn,Pn),(^n)n in [0,T] x Ss converging to 
{t,z,6). Since < Ls{zn,On) = max(a;„, L(2;„, 6'„) — e), and y„, goes to zero, this implies that 
for n large enough, x„ = Ls{zn,9n) > 0. By considering from {tn,Zn,6n) the admissible 
strategy of doing none transaction, which leads to a final liquidation value Xj^ — 
we have U{xn) < Vgitn, Zn,0n) < v{tn, Zn,9n)- Recalling Corollarv 13.11 we then obtain the 
continuity of on 5^5^ with v^{t, z,9) = U {x) = v{t, z, 9) for (z, 9) = (x, 0,p, 9) G dyS^, and 
in particular (15. Sp . Finally, we address the uniqueness issue, which is a direct consequence 
of the following comparison principle for constrained (discontinuous) viscosity solution to 

(lESD-dSSI). 

Theorem 5.2 (Comparison principle) 

Suppose u € ^-y([0,T] x Se) is a use viseosity suhsolution to (I5.5p - (l5.6p on [0, T] x Se, and 
w G ^^([0,T] X Ss) is a Isc viscosity supersolution to (|5.5p - (|5.6p on [0, T] x Se such that 

u{t,z,9) < liminf w{t' , z' ,9'), 'iit, z,9) £ [0,T] x Dq. (5.19) 

(t', z', 9') ^ {t, z, e) 
(t',z',e') e [0, T) X 5e 

T/ien, 

u < w on [0,T] X Se- (5.20) 

Notice that with respect to usual comparison principles for parabolic PDEs where we 
compare a viscosity subsolution and a viscosity supersolution from the inequalities on the 
domain and at the terminal date, we require here in addition a comparison on the boundary 
Dq due to the non smoothness of the domain 5^ on this right angle of the boundary. 
A similar feature appears also in [17], and we shall only emphasize the main arguments 
adapted from [3], for proving the comparison principle. 
Proof of Theorem 15.21 

Let u and w as in Theorem 1 5. 2 1 and (re)define w on [0, T] x dS^ by 

w{t,z,9) = liminf w{t',z',9'), {t, z,9) £ [0,T] x dSe. (5.21) 

(t', z', S') -. (t, z, 0) 
(t'.z'.S') 6 [0,T) X Se 

In order to obtain the comparison result (|5.20p . it suffices to prove that supjQ 2-]^^^(tt — w) 

< 0, and we shall argue by contradiction by assuming that 

sup (u-w) > 0. (5.22) 
[0,T]x5. 
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• Step 1. Construction of a strict viscosity super solution. 
Consider the function defined on [0, T] x Se by 

where p' > 0, and 7' G (0, 1) wih be chosen later. The function is smooth on 
[0, T) X (S/r \ Dq), and by the same calculations as in (|3.1U|) . we see that by choosing p' > 
then 

Cip > on [0,T) X {Se\Do). (5.23) 



1"7' 2^' 

dijj dip 



dt 86 
Moreover, from ()5.4p . we have 

{iP-'Hem.z^e) = eP'^'^~''^[LM{zy' -{HelMizW'] =: ^{t,z) (5.24) 

> on[0,r]x5£. 

For m > 1, we denote by 

u{t,z,d) = e^u{t,z,9), and Wm{t,z,9) = e''\w(t, z,9) + —ipit, z,9)]. 

From the viscosity subsolution property of n, we immediately see that tt is a viscosity 
subsolution to 

du du — 

mm[u - — - — - Cu , u- Heu] < 0, on [0, T) x Se (5.25) 



mm 



[u-Ul,, u-Tieu] < 0, on{r}xcS„ (5.26) 



where we set Ul^{z,9) = e^UL^{z,9). From the viscosity supersolution property of and 
the relations (I5.23p - ()5.24p . we also derive that Wm is a viscosity supersolution to 

nil) nil) 

Wm--^--^-Cwm > ^ on [0, T) X {Se \ Dq) (5.27) 

Wm-new^ > -A on[0,r]x5e. (5.28) 
m 

Wm -Ul, > on {T} X Se. (5.29) 

On the other hand, from the growth condition on u and w in ^^([0, T] x5e), and by choosing 
7' e (7, 1), we have for all (t, 9) G [0, Tf, 

lim (u — Wm)it, z,9) = —00. 

Therefore, the use function u — Wm attains its supremum on [0,T] x Se, and from (j5.22p . 
there exists m large enough, and (t, z, 9) G [0, T] x Se s.t. 

M = sup [u-Wm) = (u- Wm){i,Z,9) > 0. (5.30) 
[0,T]xSe 

• Step 2. From the boundary condition (I5.19p . we know that {z,9) cannot lie in Dq, and 
we have then two possible cases: 
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(i) iz,e) £Ss\Do 

(ii) {z,d) e dSe\Do. 

The case (i) where {z, 9) lies in is standard in the comparison principle for (noncon- 
stained) viscosity solutions, and we focus here on the case (ii), which is specific to cons- 
trained viscosity solutions. From (|5.2ip . there exists a sequence {tn, Zn, On)n>i in [0, T) x Ss 
such that 

{tn,Zn,0n,Wm{tn,Zn,6n)) > (t, Z , 6 , Wmit, Z , 9)) aS U ^ OO. 

We then set 6n = \zn — -zl + — ^1, and consider the function defined on [0, T\ x {S^f' 
by: 

^n{t,Z,d,z! ,9') = u{t,Z,e)-Wra{t,z! ,e')-Lpn{t,Z,e,z! ,9') 

ipn{t,z,9,z',9') = \t - + \z - z\^ + \9 - 91"^ 

|Z-Z'|2 + |0-0'|2 /d{z',9') \4 
+ 2fn + \d{Zn,9n) " ^) ' 

Here, d{z,9) denotes the distance from {z,9) to dSs- Since {z,9) ^ Dq, there exists an 
open neighborhood V of {z,9) satisfying V Ci Dq = ^, such that the function d{.) is twice 
continuously differentiable with bounded derivatives. This is well known (see e.g. |12j ) 
when {z, 9) lies in the smooth parts of the boundary dS^ \ (Di^e U D2^e)- This is also true 
for {z,9) £ Dk^e for ^ £ {1)2}. Indeed, at these corner lines, the inner normal vectors 
form an acute angle (positive scalar product), and thus one can extend from {z,9) the 
boundary to a smooth boundary so that the distance d is equal, locally on the neighborhood, 
to the distance to this smooth boundary. From the growth conditions on u and w in 
^^([0, T] X 5e), there exists a sequence {tn, Zn,9n, z'^-, 9'^^) attaining the maximum of the use 
on [0, T] X {SeY- By standard arguments (see e.g. [3] or [17]), we have 

{in,Zn,9n,Zn,0'n) {t,Z,9,Z,9) (5.31) 



25n ^d{Zn,i 

u{in,Zn,9n) - Wm{in,z'^,On) ' {u - Wm){t, Z , 9) . (5.33) 



The convergence in (15.32P shows in particular that for n large enough, d{zn, 9'^) > d{zn, 9n)l'l 
> 0, and so {z'^,9'^) E Sg. From the convergence in (j5.3ip . we may also assume that for 
n large enough, {zn,9n), {z'^-,Q'n) lis in the neighborhood V of {z,9) so that the derivatives 
upon order 2 oi d{.) at {zn,9n) and {z'^,9'^) exist and are bounded. 

• Step 3. We show that for n large enough, 

u{in,Zn,9n) -UeU{in,Zn) > 0. (5.34) 

Otherwise, up to a subsequence, we would have for all n: 

u{in,Zn,9n) -'HeU{in,Zn) < 0. 



32 



By sending n to infinity, and from the upper-semicontinuity of TisU, we get with (j5.3ip : — oo 
< u{t,z,9) < 7isu{t, z,9), which shows in particular that Cs{z,9) is not empty. Moreover, 
by the viscosity supersolution property (j5.28p . we have 

m 

By substracting the two previous inequahties, we would get 

'^itni Zn ) 0n) \tn i Z^, 0^) <J 'hLeUitn, Zn) "HsWjn (tn, Z^ , 9^ Z,^^ ^n) ' 

m 

By sending n to infinity, and from the upper-semicontinuity of TieU, the lower-semicontinuity 
of Ti-eWm and A, this yields with (jS.Sip . (j5.33p 

{u-wm)it,z,e) < nsuii,z,e)-neWm{t,z,e)-—A{t,z,e). 

m 

Now, by compactness of Ce{z,9) / 0, there exists e G Cs{z,9) such that 7isu{t, z,6) = 
u{t,Te{z,6,e),0) and so 

M = (u-Wm){i,z,e) < u{t,Te{z,e,e),0)-Wm{i,re{z,e,e),0) - —A{t,z,9) 

m 

< M - —A(i,z,9), 
m 

a contradiction. 

• Step 4- We check that, up to a subsequence, tn < T for all n. On the contrary, tn = i = 
T for n large enough, and we would get from ()5.34p and the viscosity subsolution property 

u{T,Zn,9n) < UhAZn^On)- 

Moreover, by (I5.29|) . we have Wm{T, z'^^O'^) > Ui^izn, 9'^^), which combined with the former 
inequality, implies 

u(T,Zn,9n) - WrniT,z'^,0'J < UL,{Zn, 9n) - U L^z'ti, ^'n) ■ 

By sending n to infinity, this yields with (|5.3ip . (j5.33p and continuity oi Ul^: Ad = [u — 
Wm)it-,z,9) < 0, a contradiction with (j5.30p . 

• Step 5. We use the viscosity subsolution property (j5.25p of u at Zn, 9n) G [0, T) x 5^, 
which is written by (|5.34p as 

dxi du 

[U - — - — - Cu){tn.,Zn,9n) < 0. (5.35) 

The above inequality is understood in the viscosity sense, and applied with the test function 
{t,z,9) ipn{t, z,9, z!^,9'^), which is in the neighborhood [0,r] x V of {in,Zn,9n)- We 
also write the viscosity supersolution property (j5.27p of Wm at (t„, z^, 9'^) £ [0, T) x {Se\Do): 

(Wm - CWm)itn,z'n,9'J > 0. (5.36) 

The above inequality is again understood in the viscosity sense, and applied with the test 
function {t,z',9') — > —ipn{t,Zn,9n,z',9'), which is in the neighborhood [0, T] x V of 
{tn, z'^,9'^). The conclusion is achieved by arguments similar to [T7|: we invoke Ishii's 
Lemma, substract the two inequalities (j5.35p - (j5.36p . and finally get the required contradic- 
tion M < by sending n to infinity with dOTll - ([5:32] ) - (103]) . □ 
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6 An approximating problem with utility penalization 



We consider in this section another perturbation of our initial optimization problem by 
adding a cost e to the utility at each trading. We then define the value function Ve on 
[0,r] x5by 



Ve{t,Z,9) 



sup E 



UL{ZT,eT) -eNria) , {t, z,9) € [0,T] x S. (6.1) 



The convergence of this approximation is immediate. 

Proposition 6.1 The sequence (v^)^ is nondecreasing and converges pointwise on [0, T]xS 
towards v as e goes to zero. 

Proof. It is clear that the sequence (v^)^ is nondecreasing and that Vs < v on [0, T] x S 
for any e > 0. Let us prove that lim^^o = v. Fix n G N* and {t, z, 9) G [0, T] x S and 
consider some a^") G A\{t., z, 9) such that 



E 



[/l(4"\0t^) > v{t,z,9)--, 



n 



where (Z^^^G^*^)) is the associated controlled process. From the monotone convergence 
theorem, we then get 



\\mve{t,z,9) > E UL[zi'\@P 



> v{t,z,9) 



1 

n 



□ 



By the arbitrariness of n G N*, we conclude that lim^Vs > v, which ends the proof since 
we already have < v. 

The nonlocal impulse operator associated to (j6.ip is given by 

neip{t,z,9) = nip{t,z,9)-e, 

and we consider the corresponding dynamic programming equation: 



dw dw 
mm - - - , w - HeW 
ot o9 

min [w — Ul , w — HgW 



0, in [0, T) X S, 
0, in {T} X S. 



(6.2) 
(6.3) 



By similar arguments as in Section 5, we can show that Vs is a constrained viscosity 
solution to (j6.2p - (j6.3p . and the following comparison principle holds: 

Suppose u G G'yiiO, T] x 5) is a use viscosity subsolution to (|6.2p - (|6.3p on [0, T] x S, and w 
G ^^([0,T] X 5) is a Isc viscosity supersolution to (j6.2p - (j6.3p on [0,T] x S, such that 

w{t',z',9'), y{t,z,9) €[0,T]x Do. 



u{t,z,9) < 



lim inf 

{t',z', 0') (t, z, 9) 
{t',z',e') 6 [0,T) X S 



Then, 



u < w on [0, T] X S. 



(6.4) 



The proof follows the same lines of arguments as in the proof of Theorem 15.21 (the function 
il) is still a strict viscosity supersolution to (j6.2p - (l6.3p on [0, T] x S), and so we omit it. 
As a consequence, we obtain a PDE characterization of the value function v. 
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Proposition 6.2 The value function v is the minimal constrained viscosity solution in 
Q^{[0,T] X S) to (|4.4p - (|4.5p . satisfying the boundary condition 

lim v(t',z',9') = v(t,z,9) = U(0), y(t, z,e) e [0,T] x Dq. (6.5) 
(t',z',e')~>(t,z,e) 



Proof. Let V G Q-y{[^,T] x cS) be a viscosity solution in Q^{[Q,T] x S) to (|i:il) - (|i3D . 
satisfying the boundary condition (|6.5|) . Since 7i > Tie, it is clear that is a viscosity 
supersolution to (|6.2|) - (|6.3|) . Moreover, since lim(^/ 2',0')^{t,z,e) ^tit' , z' ,9') = C/(0) = v(t,z,9) 
> v*{t,z,9) for {t,z,9) £ [0,T] x Dq, we deduce from the comparison principle (j6.4p that 
^ > > > on [0, r] X 5. By sending e to 0, and from the convergence result in 
Proposition l6.lt we obtain: V > v, which proves the required result. □ 

Appendix: constrained viscosity solutions to parabolic QVIs 

We consider a parabolic quasi-variational inequality in the form: 

+ F{t,x,v,D^v,Dlv) , V -Hv] = 0, in[0,r)xO, (A.l) 



mm 



dt 

together with a terminal condition 

[v-g , v-Tiv] = 0, in {T} x O. (A.2) 



mm 



Here, O C M*^ is an open domain, F is a continuous function on [0, T] x M'' x M x M*^ x 5*^ 
(5"^ is the set of positive semidefinite symmetric matrices in M'^^'^), nonincreasing in its last 
argument, g is a continuous function on O, and 7^ is a nonlocal operator defined on the set 
of locally bounded functions on [0, T] x O by: 

7{v{t,x) = sup [v{t,T{t,x,e)) + c{t,x,e)]. 

eeC(t,x) 

C{t,x) is a compact set of a metric space E, eventually empty for some values of (t, x), in 
which case we set TCv{t,x) = —0, and is continuous for the Hausdorff metric, i.e. if {tn,Xn) 
converges to {t,x) in [0, T] x O, and (e„) is a sequence in C(t„,a;„) converging to e, then 
e S C{t,x). The functions T and c are continuous, and such that T(t,x,e) £ O for all e S 
C{t, x, e). 

Given a locally bounded function u on [0,T] x O, we define its lower-semicontinuous 
(Isc in short) envelope ti* and upper-semicontinuous (use) envelope u* on [0, T]x S by: 

Utf{t,x) = liminf u(t',x'), u*{t,x) = limsup u{t',x'). 

{t',^') ^ (t,x) (t',x') -> {t,x) 

) e [o,T) X o (t',x')elo,T)xO 

One can check (see e.g. Lemma 5.1 in |17j ) that the operator 7i preserves lower and upper- 
semicontinuity: 

(i) TCu^ is Isc, and TCu^ < (TCu)^, (ii) TCu* is use, and iTLu)* < TLu* . 

We now give the definition of constrained viscosity solutions to (IA.ip - (|A.2p . This notion, 
which extends the definition of viscosity solutions of Crandall, Ishii and Lions (see [lOj). 
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was introduced in [27j for first-order equations for taking into account boundary conditions 
arising in state constraints, and used in [29] for stochastic control problems in optimal 
investment. 

Definition A.l A locally bounded function v on [0,T] x O is a constrained viscosity solu- 
tion to (|A.ip - ()A.2p if the two following properties hold: 

(i) Viscosity supersolution property on [0,T] x O: for all {t,x) G [0, T] x O, and if S 
C^'^([0,T] X O) with = (v* - ^p){t,x) = min(t;* -(/?), we have 

dip 



mm 



dt 



{i,x) + F{i,x,(f^{i,x),Dri:(p{i,x),D^(f{i,x)) , 

v^(t,x) -nv,(t,x) > 0, (t,x)£[0,T)xO, 
[v4t,x)-g{x) , v4t,x) -nv^{t,x)\ > 0, {t,x) G {T} x O. 



mm 



(a) Viscosity subsolution property on [0, T] x O: for all {t,x) £ [0, T] x O, and ip S 
C^'2([0,T] X O) with = {v* - ip){i,x) = max(z;* -ip), we have 



mm 



dip 
'dt 



+ F{t, X, ip^{t, x), D^ip{i, x), Dlip{t, x)) , 

v*{t,x) -nv*{i,x)\ < 0, {i,x) e[0,T) xO, 
mm[v*{t,x) - g{x) , v*{t,x) -nv*{t,x)\ < 0, {t,x)£{T}xO. 
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